

AERONAUTIC SYMBOLS 

1. FUNDAMENTAL AND DERIVED UNITS 


Metric 


English 


Symbol 


Length, 
Time,_- 
Force , ... 


meter 

seconds _ 

weight of 1 kilogram. 


Power 

Speed 


Y I kilometers per hour. 

4 \ meters per second _ _ 


Abbrevia- 

tion 

Unit 

Abbrevia- 

tion 

I 

m 

i 

foot (or mile) 

ft (or mi) 

s 

! second (or hour) 

sec (or hr) 

kg 

weight of 1 pound 

lb 

! --! 

horsepower _ _ ... 

hp 

• kph i 

[ miles per hour ! 

j mph 

j *nps 

feet per second - 

j 


Weight —mg 

Standard acceleration of gravity— 9. 80665 m/s 2 
or 32.1740 ft/sec 2 
W 

Mass— — - 

g 

Moment of inertia --m/r. (Indicate axis of 
radius of gyration k by proper subscript.) 
Coefficient of viscosity 


2. GENERAL SYMBOLS 

> Ki 


v Kinematic viscosity 

p Density (mass per unit volume) 

Standard density of dry air, 0.12497 kg-m“ 4 -s 2 at 15° C 
and 760 mm; or 0.002378 ih-ft' 4 sec 2 
Specific weight of “standard" air, 1.2255 kg/rn 3 or 
0.07651 Ib/ou ft 


3. AERODYNAMIC SYMBOLS 


s 

Area 

iw 

Angle of setting of wings (relative to thrust line) 

Su 

Area of wing 

it 

Angie of stabilizer setting (relative to thrust 

G 

Gap 


line) 

b 

Span 

Q 

Resultant moment. 

c 

Chord 

11 

Resultant angular velocity 

A 

b 2 

Aspect ratio, ^ 

R 

VI 

Reynolds number, p — where i is a linear dimcn- 
- y 

V 

True air speed 
Dynamic pressure, 


sion (e.g., for an airfoil of 1.0 ft chord, 100 mph, 
standard pressure at 15° C. the corresponding 

V 


Reynolds number is 935,400; or for an airfoil 

L 

j. 

Lift, absolute coefficient 

’ qS 


of 1.0 m chord, 100 mps, the corresponding 
Reynolds number is 6,865,000) 

D 

D 

Drag, absolute coefficient C D —- r , 

a 

Angle of attack. 


qS , 

€ 

Angle of downwash 

Do 

Profile drag, absolute coefficient C D 

a 0 

Angle of attack, infinite aspect ratio 



Angle of attack. Induced 

D x 

f) 

Induced drag, absolute coefficient C Di ~—f 

( A' 

Kg 

Angle of attack, absolute (measured from zero- 
lift position) 

A 

a 

Parasite drag, absolute coefficient C Dp — -J 

G 

Cross-wind force, absolute coefficient 

y 

Flight-path angle 



N 
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GENERAL THEORY OF AERODYNAMIC INSTABILITY AND THE MECHANISM OF 

FLUTTER 

By Theodore Theodorsen 


SUMMARY 

The aerodynamic forces on an oscillating airfoil or 
airfoil-aileron combination of three independent degrees 
of freedom have been determined. The problem resolves 
itself into the solution of certain definite integrals, which 
have been identified as Bessel functions of the first and 
second kind and of zero and first order. The theory, 
being based on potential flow and the Kutta condition, 
is fundamentally equivalent to the conventional wing- 
section theory relating to the steady case. 

The airforces being known, the mechanism of aerody- 
namic instability has been analyzed in detail. An exact 
solution, involving potential flow and the adoption of the 
Kutta condition, has been arrived at. The solution is of 
a simple form and is expressed by means of an auxiliary 
parameter k. The mathematical treatment also provides 
a convenient cyclic arrangement permitting a uniform 
treatment of all subcases of two degrees of freedom. The 
flutter velocity, defined as the air velocity at which flutter 
starts , and which is treated as the unknown quantity, is 
determined as a function of a certain ratio of the fre- 
quencies in the separate degrees of freedom for any magni- 
tudes and combinations of the airfoil-aileron parameters. 

For those interested solely or particularly in the numeri- 
cal solutions Appendix I has been prepared. The rou- 
tine procedure in solving numerical examples is put 
down detached from the theoretical background of the 
paper. It first is necessary to determine a certain number 
of constants pertaining to the case, then to perform a few 
routine calculations as indicated. The result is readily 
obtained in the form of a plot of flutter velocity against 
frequency for any values of the other parameters chosen. 
The numerical work of calculating the constants is sim- 
plified by referring to a number of tables , which are in- 
cluded in Appendix I. A number of illustrative examples 
and experimental results are given in Appendix II. 

INTRODUCTION 

It has been known that a wing or wing-aileron struc- 
turally restrained to a certain position of equilibrium 
becomes unstable under certain conditions. At least 
two degrees of freedom are required to create a con- 
dition of instability, as it can be shown that vibrations 


of a single degree of freedom would be damped out by 
the air forces. The air forces, defined as the forces due 
to the air pressure acting on the wing or wing-aileron 
in an arbitrary oscillatory motion of several degrees of 
freedom, are in this paper treated on the basis of the 
theory of nonstationary potential flow. A wing- 
section theory and, by analogy, a wing theory shall be 
thus developed that applies to the case of oscillatory 
motion, not only of the wing as a whole but also to 
that of an aileron. It is of considerable importance 
that large oscillations may be neglected; in fact, only 
infinitely small oscillations about the position of 
equilibrium need be considered. Large oscillations 
are of no interest since the sole attempt is to specify 
one or more conditions of instability. Indeed, no 
particular type or shape of airfoil shall be of concern, 
the treatment being restricted to primary effects. The 
differential equations for the several degrees of freedom 
will be put down. Each of these equations contains a 
statement regarding the equilibrium of a system of 
forces. The forces are of three kinds: (1) The inertia 
forces, (2) the restraining forces, and (3) the air forces. 

There is presumably no necessity of solving a general 
case of damped or divergent motion, but only the 
border case of a pure sinusoidal motion, applying to the 
case of unstable equilibrium. This restriction is par- 
ticularly important as the expressions for the air force 
developed for oscillatory motion can thus be employed. 
Imagine a case that is unstable to a very slight degree; 
the amplitudes will then increase, very slowly and the 
expressions developed for the air forces will be appli- 
cable. It is of interest simply to know under what 
circumstances this condition may obtain and cases in 
which the amplitudes are decreasing or increasing at a 
finite rate need not be treated or specified. Although 
it is possible to treat the latter cases, they are of no 
| concern m the present problem. Nor is the interna! 
or solid friction of the structure of primary concern. 
The fortunate situation exists that the effect of the 
solid friction is favorable. Knowledge is desired con- 
cerning the condition as existing in the absence of the 
internal friction, as this case constitutes a sort of lower 
j limit, which it is not always desirable to exceed. 
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Owing to the rather extensive field covered in the 
paper it lias been considered necessary to omit many 
elementary proofs, it being left to the reader to verify 
certain specific statements. In the first part of the 
paper, the velocity potentials due to the flow around 
the airfoil-aileron are developed. These potentials 
are treated in two classes: The noncirculating flow 
potentials, and those due to the surface of discon- 
tinuity behind the wing, referred to as “ circulatory ” 
potentials. The magnitude of the circulation for an 
oscillating wing-aileron is determined next. The 



Figure 1.— Conformal representation of the wing profile by a circle. 

forces and moments acting on the airfoil are then 
obtained by integration. In the latter part of the 
paper the differential equations of motion are put 
down and the particular and important case of un- 
stable equilibrium is treated in detail. The solution 
of the problem of determining the flutter speed is 
finally given in the form of an equation expressing a 
relationship between the various parameters. The 
three subcases of two degrees of freedom are treated 
in detail. 

The paper proposes to disclose the basic nature of 
the mechanism of flutter, leaving modifications of the 
primary results by secondary effects for future investi- 
gations. 1 Such secondary effects are: The effects of a 
finite span, of section shape, of deviations from poten- 
tial flow, including also modifications of results to 
include twisting and bending of actual wing sections 
instead of pure torsion and deflection as considered in 
this paper. 

The supplementary experimental work included in 
Appendix II similarly refers to well-defined elementary 
cases, the wing employed being of a large aspect ratio, 
nondeformable, and given definite degrees of freedom 
by a supporting mechanism, with external springs 
maintaining the equilibrium positions of wing or wing- 
aileron. The experimental work was carried on 
largely to verify the general shape of and the approxi- 
mate magnitudes involved in the theoretically pre- 
dicted response characteristics. As the present report 
is limited to the mathematical aspects of the flutter 
problem, specific recommendations in regard to prac- 
tical applications are not given in this paper. 

1 The effect of internal friction is in some cases essential; this subject will be 
contained in a subsequent paper. 


VELOCITY POTENTIALS, FORCES, AND MOMENTS OF 
THE NONCIRCULATOR Y FLOW 


We shall proceed to calculate the various velocity 
potentials due to position and velocity of the individ- 
ual parts in the whole of the wing-ailerOQ system. 
Let us temporarily represent the wing by a circle (fig. 
1). The potential of a source <= at the origin is given 

by 

l°g (tf + y*) 

For a source c at (x h y t ) on the circle 


<p= ir log {( x ~ x i) 2+ (y~yty} 

Putting a double source 2c at (xi,y\) and a double 
negative source —2c at (xi,~ 2 /i) we obtain for the flow 
around the circle 



(g~gQ a + ( y-yi Y 

(x-zi ) 2 +(y+yi ) 2 


The function <p on the circle gives directly the sur- 
face potential of a straight line pq, the projection of the 
circle on the horizontal diameter. (See fig. 1.) In 


this case y = ^\ — x 2 and cp is a function of x only. 
We shall need the integrals: 



(x-xr ) 2 +(y-yQ 2 
(x-xy+iy+yj 2 


dx 1 ==2(x-c)\ogN—2-y/l—x i cos x c 


and 

/. 


l[og (»- xoffiy + $ - vr^vi 3 ? 


where 


■cos X c{x— 2c) -yjl — 2 ^+ (x — c) 2 log N 
1 ~Z?tJ 1 —C 2 


The location of the center of gravity of the wing- 
aileron x a is measured from a, the coordinate of the 
axis of rotation (fig. 2); x& the location of the center 



Figure 2. — Parameters of the alrfoil-aUeron combination. 


of gravity of the aileron is measured from c, the coordi- 
nate of the hinge ; and r a and r$ are the radii of gyration 
of the wing-aileron referred to a , and of the aileron 
referred to the hinge. The quantities x$ and ra are 
‘‘reduced ;; values, as defined later in the paper. The 
quantities a, x a , c, and are positive toward the rear 
(right), h is the vertical coordinate of the axis of rota- 
! tion at a with respect to a fixed reference frame and is 
positive downward. The angles a and j3 are positive 
clockwise (right-hand turn). The wind velocity v is to 
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the right and horizontal. The angle (of attack) a 
refers to the direction of v, the aileron angle refers to 
the undeflected position and not to the wind direction. 
The quantities r a and r 0 always occur as squares. 
Observe that the leading edge is located at — 1, the 
trailing edge at +1. The quantities a , c, x a} x 0 , r a , 
and />, which are repeatedly used in the following 
treatment, are all dimensionless with the half chord b 
as reference unit. 

The effect of a flap bent down at an angle /S (see fig. 
2) is seen to give rise to a function <p obtained by sub- 
stituting — vpb for «; hence 

<Pq = —jjr [ -yfl -X 2 COS^'C- ( X—C ) logiVj 

To obtain the effect of the flap going down at an 
angular velocity 0, we put e = — (x x — c)/3& 2 and get 

n = cVl - * 2 + c °s-*c(* - 2c) V 1 - r* 

— (x — c) 2 log Af] 

To obtain the effect of an angle a of the entire air- 
foil, we put c= — 1 in the expression for <p 0 , hence 
<p a — votb Vl — x 2 

To depict the airfoil in downward motion with a veloc- 
ity h (4- down), we need only introduce ~ instead of a. 

Thus _ 

(p^ — hb-yjl—x 2 

Finally, to describe a rotation around point a at an 
angular velocity a, we notice that this motion may be 
taken to consist of a rotation around the leading edge 
c— — 1 at an angular velocity a plus a vertical motion 
with a velocity — a(l +a)b. Then 

<Pa — — 7r(x4- 2) Vi — x 2 — «(1 4 a)6 2 Vl ~ x 2 

Ztt 

= ctb 2 (j^x — a ^ Vi ” x 2 


INTEGRALS 


j <P« dx 

J c 

X 

jc 

J* (j P 0 dx 
<Pa(x- 

jc <Ph ^' X ~ 

fe 

jc <P ^ X ~ 

Jc <P ^ X ~ 


- % VaT< 
- \hT \ 


= ab 2 Ta 




V 2 


pt 2 


c)dx = — -^oaT x 


c)dz = — ^ hT x 


c)dx = cib 2 T n 


c)dx = --j-v$T 2 


c)dx=-~jT 3 
47 r 


f+> , b 

I ip a dx — ^vair 
J x (phdx^^T 


. ,,7ra 
ab -g 


\v»T< 


X+l 

^ <t pidx = 

f+i 

J ^ <psdx = 

J t •Pbdx— - ~PT X 

J'_ i <P*(x-c ) dx = - 
J* <ph(x- 

r+i 

I <pz(x — c)dx = ab 2 TuTr 

C+i b 

I <Pti(x-c)dx « ~2 vQTf 


^Vacw 


c) dx = — ^ Ac 7 t 


CONSTANTS 

T\ = — gVl — C 2 (2 4 c 2 ) +C COS^C 

T 2 = c(l — c 2 ) — VW(1 4c 2 )cos _1 c4-c(cos~ ! c) 2 

T 3 = — ^- + c 2 ^ (cos _1 c) 2 + “C Vl cos _1 c(7 4 2c 2 ) 

- |(l-c 2 ) (5c 2 + 4) 

Tt = — cos _1 c 4 c Vl — c 2 

T 5 = — (1 — c 2 ) — (cos -1 c) 2 4 2c Vl — c 2 cos~‘c 

r 6 =r 2 

T 7 = — (^| + c2 ^) COS _1 C + ^C -yJT—C 2 (7 + 2c 2 ) 


The following tables give in succession the velocity 
potentials and a set of integrals 2 with associated con- 
stants, which we will need in the calculation of the air 
forces and moments. 

VELOCITY POTENTIALS 
ip a — Vab V 1 ~X 2 
tpk^hb^l — x 2 

<Pi = a& 2 Q z - a) V 1 “ sc 2 

<P 0 = --i>/S&[Vl — z 2 cos - 1 c-~ (x — c) log A 7 ] 

t r 

= OTVi — ri*Vi — r ! + (x — 2c) \T — ri'-cos~ I c 

47T 

— (x — c ) 2 log N) 

where 1 - C x- Vi^^Vl- c 2 

x — c 

* Some of the more difficult integral evaluations are given ki Appendix III. 


Tg == — )^^l—c 2 (2 c 2 4- 1) 4- c cos -1 c 

r 9 = ^[^(Vw) + ^]==^(-2) + a2 T 4 } 

where ^ | ^ Vl — c 2 ^) 3 

T w = Vl ~-c 2 4-cos _1 c 

Tjj = cos -1 c (1 — 2c) 4 Vl — c 2 (2 — c) 

Tn— Vl"C 2 (2 4c)~ cos -1 c (2c 4 1) 


Tu = 7 ; ( — T 7 — (c — a) T ,1 


r, , a= -- 


The velocity potentials being known, we are able to 
calculate local pressures and by integration to obtain 
the forces and moments acting on the airfoil and 
aileron. 
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Employing the extended Bernoulli Theorem for un- 
steady flow, the local pressure is, except for a constant 


(w 2 '0<p\ 

Pf-e^ + dt) 


where w is the local velocity and <p the velocity poten- 
tial at the point. Substituting = + ^ we obtain 

ultimately for the pressure difference between the 
upper and lower surface at x 

P = ~ 2p ( v Tx + W 

where v is the constant velocity of the fluid relative to 
the airfoil at infinity. Putting down the integrals for 
the force on the entire airfoil, the moment on the flap 



~=5~ 

-AT 


' Figure 3. — Conformal representation of the wing profile with reference to the 
circulatory flow. 

around the hinge, and the moment on the entire air- 
foil, we obtain by means of partial integrations 

J +i 

tpdx 

Mp = — 2pb 2 j <p(x — c) dx + 2pvb ( <pdx 
Jc Je 

r + 1 p+i 

M a = —2pb 2 I <p(x — c)dx + 2pr&J ^ <pdx 

r + 1 

— 2p6 2 I <p{c — a)dx 

Or, on introducing the individual velocity potentials 
from page 5, 

P— — pb 2 [rjra+ Trh — bra-ix — vT $ — bT\$\ (I) 

Mp~ — pb 2 — vT\a— Tih+ 2l\zba — ^ ^ T 3 bfi J 

+ pvb 2 j — vT t a — Tji + 2T< j ba — — vJ\ ft — ^ 7 \bji J 
= -pb^T^a- (2T,+ Ti)bva + 2T n b 2 d+ ^ 7> 2 0 
+ (~ T 2 -± bvfl - i b 2 Tj + Tffih - I'M l J (II) 

pb 2 £ 7r } '~ (x i 4- tr ) h 2 c i f v : TS \ 1 1\ - T 

- (c — a) } bfio + { — T 7 — (c - a) 2 \ } 6 2 £ 

— buiri , — Trr/t, J (III) 




VELOCITY POTENTIALS, FORCES, AND MOMENTS 
OF THE CIRCULATORY FLOW 

In the following we shall determine the velocity 
potentials and associated forces and moments due to a 
surface of discontinuity of strength U extending along 
the positive x axis from the wing to infinity. The 
velocity potential of the flow around the circle (fig. 3) 
resulting from the vortex element — AF at ( X 0 , 0) is 

Y 


af 

<f>r ~ 27T 


tan' 


X-Xo 


— tan 


, Y I 

*-v 


Ar * -1 

“2i ta " 


(-Vo + ^) : 


X 2 




X+Y 2 + 1 


where (A r , Y) are the coordinates of the variable 
and X 0 is the coordinate of — Ar on the x axis. 

1 

Introducing X 0 + y- — 2x 0 
Ao 

or Ao=x 0 + a/xo 2 — 1 on the x axis 
and X=x and Y=-yJl—x 2 on the circle 
the equation becomes 


AT Vl ~~ x? yjx o 2 — 1 

This expression gives the clockwise circulation 
around the airfoil due to the element — AF at x 0 . 

We hftvc:p=— 2p(^ + f^) 

But, since the element — A1 1 will now be regarded as 
moving to the right relative to the airfoil with a 
velocity v 

dtp _ c )<p 

dt dx 0 


Hence, ‘p 
Further 


) J'bv 1 dtp\ 
' p \ dx dx 0 / 


XoVl ~~ P 


2ir dtp /— a — r(l — XX 0 ) Vl — x 2 (1 — xx 0 ) 2 

AFdx“ Vr ° 1 


1 + 


V*o a 


(1 -x 2 )(xq 2 - 1) 

(1 - xx 0 y 

1 


and 


Vl-x a (*o “ x) 

1 X 0 


2 7T dtp 




i(l-XXo) \x 0 2 - 1 . 


+ V*o 2 - 1 


(1 — XXq) 2 


AF dx 0 


By addition: 


1 + 


(1 x 2 ) (Xq 2 1) 


(1 XX 0 ) 2 

\ r :r F 1 

v j: 0 - 1 0 3") 


dtp t dtp AT 
t — 


x 0 +x 


dx dx 0 2k 71 — P-yjx^ — 1 
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■/ 


To obtain the force on the aileron, we need the 
integral 

f Y ** + dx - — f 2a±* dx 

X \dx + dr J ax 2 X Jc Vx„ J - 1 VF-T= 01 


Arf j 0 




2tLV^o 2- ' 1 


V 1 — s' 
Vzo 2 - 


.2 “ll 

COS" l X + 

1 Jc 


Arf j 0 , Vi — c 1 

= 2iLvW 0S c+ vi 


— 1 

0 J -1J 


Thus, for the force on the aileron 




AP c i = - pvb 


AP,i = 


=cos ! cH — t 
1 -y/^o' 




VVV 


py&“ — £4 == (cOS j1 C — V~1 "— C 2 ) 

T LV¥~ 1 


or 


+ Vjfi V i ” c '] 

Integrated, with AT = C7dj 0 

F el = -^[(cos-c- Vl — c 2 )J^ 

+ Vl -~ 


xo 


V^O 2 - 1 


= £/dx 0 


for c~- — 1 we obtain the expression for P , the force 
on the whole airfoil 


— r 


X 0 


Vxo 2 


C/d j 0 


(IV) 


Since U is considered stationary with respect to the 
fluid elements 

U=l(vt — J 0 ) 

where f is the time since the beginning of the motion. 
U is thus a function of the distance from the location 
of the first vortex element or, referred to a system 
moving with the fluid, U is stationary in value. 

Similarly we obtain for the moment on the aileron 


sx 


dip j d<p 

dx dj 0 , 

Ar 


\ / X , Ar P (x-c) (j 0 + j) , 
J(z-e)d^^J c 


2t Vxo 1 


FiD° Vi 


T— — 2 , xVl-X 2 


J 2 + 


cVl -X 2 


+ Q- J^cos-j] 

= + ^VS7^[ (Xo+ 2" c)Vi ~ 

+ — 2joC^COS _1 cJ 


+ 


A! 


.. .if Y r 

LvV-i' 

l 


t) 


+ 






cos” 'c — r 


Vl-c*) 


Fincfly 

A^=-M^[^ rT {Vw(l + |) 

cos -1 c(c +|)} +5-y/^rl (cos-'c-cV'F 1 ?) J 
Putting Ar= C7dj 0 and integrating 

. M *— fTP^O+C) 


— cos c 


H)ii: 


x 0 


Udx 0 


yxo 2 - 1 

+ (cos _I c- cV(l ~c 2 ) 2 Ji "V^jo— 1 ^ 

Further, for the moment on the entire airfoil around a 

Hi&sO ( - t ” a)dx= ” lx viki[( x " + l ” “) vrr? 
+ (i-^) cos ”' x T-'r 

1 

2 — Xo& 

and AM, = — pvb 2 AT J^Tz Tf 

Integrated, this becomes 


pm J 

„ I 2 
JS= t' 

nth 

Ji . Vx 0 2 ~ 1 


= — pvb 2 
= — pvb 2 


< tt+ £) | 

Vj 0 2 -1 I 


C/dj 0 



T 

-I 12' 


X 0 + 


T~( a + I 


\ x c 


c/dj„ (VI) 


Jn-I \ ^/Vjo 2 -1) 

THE MAGNITUDE OF THE CIRCULATION 


The magnitude of the circulation is determined by 
the Kutta condition, which requires that no infinite 
velocities exist at the trailing edge, 
or, at x— 1 

~(<p r + <p a + <pk + <f>a + <pf) + <pp) = finite 
Introducing the values of <p a , etc. from page 5 and 

N * 

<P T from ^ page 6 gives the important relation: 


1 r m Vxo+i 
Vx o~ 1 


2irji 


Udx 0 = va + h + b 


(i-a)« 


+ —V0 + 6—4 

nr '.nr 


(VII) 


This r - 5 at ion tmmt !: ,m / to comply with t he 

Kutta condition, v. .uvi) . ■ > .inti the flow shall leave 

the airfoil at the trailing edge. 

It is observed that the relation reduces to that of the 
Kutta condition for stationary flow on putting j„= oo, 
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and in subsequence omitting the variable parameters 
d, $, and 
Let us write 


Udxo=v a + h+b(±-i)« 
Tn 




Introduced in (IY) 

P — — 2vpvbQ 

from (V) 








— 2pvb 2 


(^Jl-c 2 (l+^~cos 1 c(c + £)) X 


nm^- ( ’ 


Q 



and from (VI) 
M a — — 2tt pvb 2 
Introducing 


1 / , l\Ji V^o 2 “l 
2“( a+ 2 


Udx 0 


4 _ 


Q 


s: 


x 0 


W- 1 


Udxo 


Udx 0 


we obtain finally 


where Q is given above and C~ C(k ) will be treated in 
the following section. 

VALUE OF THE FUNCTION C (*) 

Put 

where s = vt (s— »«>), the distance from th z first vortex 
element to the airfoil, and k a positive constant deter- 
mining the wave length, 
then 


r fa> +i 

Ji v* 0 -i 

P=~2pvbicCQ (VIII) 

— 2 pvb 2 ^ Vl ~ c 2 ^l 4- -cos -1 c(c 

+ ~(cos -1 c-cVW)]^ -pvb 2 {T l2 C-l\)Q (IX) 
M a = 2*pvb 2 [(a + ^ C - ~] (2 (X) 


<**)■ 


j; 


x 0 


V^O 2 ' 


e ikx odxn 


s; 


x ° ^ \e ikx odxa 


(XI) 


V*o 2 ~ 1 

These integrals are known, see next part, formulas 
(XIV) — (XVII) and we obtain 3 


C{k) = 




- J x + iY x 


T r 7T„ , -TTtz T -(Jl + PoHiOWo) 

(- Ji+iFi)[ - ( Ji +Y 0 )-i(Y l - J 0 )} 
(J l +Fo) 2 +(Vi-</o ) 2 

t/i(Ji+r 0 ) + FiOWo) 

' (J 1 +Fo) 2 +(V 1 -Jo ) 2 

. ^(J^r^-J^n-Jo) _ p . r 
' l (j l +Y 0 y+(Y l -j 0 y r+l{J 


where 


F= 

G~ 


J\(Ji + Fq) + Fi(Fi~ Jq) 

(Ji + Fo) 2 +(Fi-»/o ) 2 

YrYo + J.Jo 


(XII) 

(XIII) 


(J l +F 0 ) 2 +(F l -Jo ) 2 
These functions, which are of fundamental import- 
ance in the theory of the oscillating airfoil are given 

graphically against the argument 4 m figure 4. 

SOLUTION OF THE DEFINITE INTEGRALS IN C BY MEANS OF BESSEL 
FUNCTIONS 

We have 


i&w-jr 


g-* C 08 h< pQgk n t ^ 


where 

and 

but 


(Formula (34), p. 51— Gray, Mathews 
& MacRobert: Treatise on Bessel 
Functions. London, 1922) 

inr 

K n (t)=e 2 G„ (it) 

(Eq. (28), sec. 3, p. 23, same reference) 
<7, {x)--Y, (x) + [log 2 - T + ^]j„ (i) 


Yn (Z) “ | Yn (x) 4- (log 2 - 7 ) Jn ix) 


(where Y n (r) is from N. Nielsen: 
Handbuch der Theorie der Cylinder- 
funktionen. Leipzig, 1904). 

Hi(*) 


* This may also bo expressed in nankel functions, C— 


/A(»)+i Ho(») 
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Thus, 


We have 




or 


r ® c ® ptkz 

K q ( - zA') = I OOBh * d f « - r dx 
Ji Ji yx 2 — 1 

, n | . x T /, N f 09 cos &xdx , . f 

•gr.W + tgJiW-^ + 


Thus, 


I 


Further, 


I 


“ cos fcxdx 7T ... 

■^r = -2 r » (i) 

“ sin fcxdx £ 

V^ 2 — 1 2 


«/o (&) 


® sin lex dx 

V^ 2 - 1 
XXIV) 
(XV) 


Ki( — ik)=^ f e a 008,11 coshfdf = f :- X( fe 

Ji yx 2 — 1 


Thus, 


r 

r 

r 


Vx*- i 


(cos &x 4- i sin &x) dx 


x cos kxdx i r/n 

-W^r~~2 Mk) 

' x sin&xdx t 

■^r = ~2 r 'W 


(XVI) 

(XVII) 


TOTAL AERODYNAMIC FORCES AND MOMENTS 

TOTAL FORCE 

From equations (I) and (VIII) we obtain 
P— — pb 2 {vra + 7 rh~ ir bad — vT \0 — T\b0) 

— 2tp^6(7| voc 4- h 4- 6^2 — (i^a. 4- — T w v$ 

+ 6“7'n/j} (XVIII) 

TOTAL MOMENTS 

From equations (II) and (IX) we obtain similarly 
Mf = - { - 27V - 7, + T,(a - i) J vba + 2T l3 b 2 a 

+ V/3(T S - T,T k ) - ±vbpT,T n - i TM 
7)A/fj pvb 2 T iiO l^va. + 4+ fc(^ — rt)a 
+~T l0 vH+b~T n i)j (XIX) 


From equations (III) and (X) 

^ vba 4- irb 2 \ 

+ (T< 4- T l0 )v 2 p 

+ (r,-r,- (c-a)T, + ^T,^ 

~ (l) i - (c n) T i jb'fi — axi/fl 

■t- 2pvl 2 n(a + t)c j va ' ^ + a / 

+ lr lo v0+b~T u l3j (XX) 

DIFFERENTIAL EQUATIONS OF MOTION 

Expressing the equilibrium of the moments about a 
of the entire airfoil, of the moments on the aileron 
about c, and of the vertical forces, we obtain, respec- 
tively, the following three equations: 
a : — I a d — 1/0 — b(c — a)S/0 — S a h — aC a + M a — 0 

(3 • 1/0 I fid b(c ci)dSfi hSfi — fiCfi 4~ Mfi — 0 

h i hAl dS a ~ 0Sfi — hCh 4" P — 0 
Rearranged: 

a\ dl a + 0(lfi+ b(c — a)Sfi ) 4- hS a + aC a — M a = 0 

j3 : 4* b(c — &)Sfi) 4- 01 fi 4- hS$ 4- (3Cfi — Mfi = 0 

h: dS a 4- 0Sfi 4- hM 4 - KC h — P = o 

The constants are defined as follows: 


p, 

mass of air per unit of volume. 

b, 

half chord of wing. 

M, 

mass of wing per unit of length. 

S a} Sfif 

static moments of wing (in slugs-feet) per 
unit length of wing-aileron and aileron, 
respectively. The former is referred to 
the axis a; the latter, to the hinge c. 

Ictyl&l 

moments of inertia per unit length of 
wing-aileron and aileron about a and c, 
respectively. 

Cay 

torsional stiffness of wing around a, cor- 
responding to unit length. 

Cfi, 

torsional stiffness of aileron around c, cor- 
responding to unit length. 

Ch , 

stiffness of wing in deflection, correspond- 
ing to unit length. 

DEFINITION OF PARAMETERS USED IN EQUATIONS 

irpb 2 

K ~~W' 

the ratio of the mass of a cylinder of air of 
a diameter equal to the chord of the 
wing to the mass of the wing, both taken 
for equal length along span. 
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r “ _ V Mb 2 ’ 
_ S a 
Xa Mb’ 



the radius of gyration divided by b. 

the center of gravity distance of the wing 
from a, divided by b. 

the frequency of torsional vibration 
around a. 

reduced radius of gyration of aileron 

divided by b, that is, the radius at 
which the entire mass of the airfoil 
would have to be concentrated to give 
the moment of inertia of the aileron Ip. 


Xp — ■ ’ reduced center of gravity distance from c. 

!(T 

— t frequency of torsional vibration of aileron 
V J- 

around c. , . - . • 

CO/, — frequency of wing in deflection. 

FINAL EQUATIONS IN NONDIMENSION AL FORM 

On introducing the quantities M ay Mg, and P, 
replacing I\ and T l3 from page 5, and reducing to 
nondimensional form, we obtain the following system 
of equations: 


(A) aj% a ! + «(! + a 2 )J+d!*(j-a^+a-^+ii£r/ + (c-a).r fi -L 

(B) a[r,‘ + (e.— a)x t - - (c - a) A] + « (p - T, fi (r„* - 


-(c- a )^] + ^|[-2p-Q-a) 7’ 4 ] 


+ 0 


1 T 


, T u 

b + XT K 


(C) 


at Xa 


) + Aj « + $(xi, - b,« ) - $ t < k ; + h (i + «) 5 + h ji 


vC(k)[va^h + 


b L b 


G-°) 


r n 

i I U ^ j J II • 

CC+—T0 + 7,-/3 


Tip 0 

V b 1 


h 1 

Mb 


+ 2 K 


vC{k) n 
6 


Dot h 

LT + 5+ 


(r“) A+ T^ + f^]- 0 


SOLUTION OF EQUATIONS 


As mentioned in the introduction, we shall only have 
to specify the conditions under which an unstable 
equilibrium may exist, no general solution being 
needed. We shall therefore introduce the variables at 
once as sine functions of the distance s or, in complex 

form with \ as an auxiliary parameter, giving the 


ratio of the wave length to 27r times the half chord b: 


and 






a = a n e 




0=/v‘hJ + O 

A=v‘(*t + 0 


where s is the distance from the airfoil to the first 
ds 

vortex element, jy — v, and ^ and <p are phase angles 
of 0 and h with respect to a. 


Having introduced these quantities in our system of 
equations, we shall divide through by k. 

We observe that the velocity v is then contained in 
only one term of each equation. We shall consider 
this term containing v as the unknown parameter ffiY. / 
To distinguish terms containing X we shall employ a 
bar; terms without bars do not contain X. 

We shall resort to the following notation, taking care 
to retain a perfectly cyclic arrangement. Let the 
letter A refer to the coefficients in the first equation 
not containing C(k) or X, B to similar coefficients 
of the second equation, and C to those in the third 
equation. Let the first subscript a refer to the first 
variable a, the subscript 0 to the second, and h to the 
third. Let the second subscripts 1, 2, .3 refer to the 
second derivative, the first derivative, and the argu- 
ment of each variable, respectively. A al thus refers 
to the coefficient in the first equation associated with 
the second derivative of a and not containing Cik) or 
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X; Cha to the constant in the third equation attached to j 
A, etc. These coefficients 4 are as follows: ! 



A a3 — 0 


K 7T \ K 7 T / 

2 ":C -') r -] 


A s ,-\(T t +T m ) 

A.=--a 

K 

A b 2~0 



K IT ' \K T / 


B a3 = Q 


(«4«) 



The solution of the instability problem as contained 
in the system ol three equations A, B, and C is given 
by the vanishing- of a third-order determinant of com- 
plex numbers representing the coefficients. The solu- 
tion of particular subcases of two degrees of freedom 
is given by the minors involving the particular co- 
efficients. We shall denote the case torsion-aileron 
{a, /3) as case 3 , < viler on-deflectio n (ffi h) as case 2, and 
deflection-torsion (k, «) as case 1. The determinant 
form of the solution is given in the major case and in 
the three possible subcases, respectively, by: 


and 


D= 


Rtia~\ kTfraj Rq@ T' zd a @, Bah | 11 %l ah 
B bai Hbftflil bjif R (jhfl-tl frh 

Rtuflil cay Be pflilcP, Rchflilch 


Xf C h 
~Al a a~ 
M b $ = 


-0 


Baa T Haa) R<i{i~\~id a$ \ 
Rhaflil ba, Rbpflld be\ 

= 0 

Case 3 

Rb&~\~H BbhflUbh 
R cpflid c $) RcJiflil ch\ 

-0 

Case 2 

Bch~\~id-ch} Rca~\~idca\ 

B ah~\~id ah/ di aa -\-'lI aa i 

=0 

Case 1 


D 1 

£>81 =» 

K 7 T 


Pi 3^3 


Bffi 


B^^(T 5 -TJ\ o ) 

Bh-Z-J-T, 

K X 


B/i2 ~ 0 
Bh 3 = 0 

C al = ~ — o (=A„) 

O a 2=1 

O a 3 — 0 


c « = 7 ~l Tt (= ' e *' ) 

c k =--t, 

T 

Cfo~0 

C M -i + l 

K 

0,2 = 0 

0,2 — 0 

‘ The factor ^ or ^ is not included in these constants. See the expressions for 

the R's and I’a oa next page. 


REAL EQUATIONS IMAGINARY EQUATIONS 


RaaRai 3 


daad a(3 

= 0 


+ 

d aa d a8 

RbaRblS 


d bad b& 

1 ba^ 

RbaRbt 

BbpBbh 

RcpRch 

- 

db&dbh 
d c pd C h 

= 0 

Rb$Rt>h 

Icfdlch 

+ 

dbpdbh 

RcpRch 

RchRca 


d chd ca 

-0 

RchRca 

+ 

d chd ca 

B at>R act 


I ahd aa 

act 

R ahB aa 


"0 Case 3 


— 0 Case 2 
=0 Case 1 


Note. —Terms with bars contain X\ terms without bars do not contain A'. 

The 9 quantities R aa , R afi , etc., refer to the real parts 
and the 9 quantities l aa , I a g, etc., to the imaginary 
parts of the coefficients of the 3 variables a, (3, and h 
in the 3 equations A, B, C on page 10. Denoting the 
coefficients of a, a, and « in the first equation by p , 
q, and r, 


' B aa + id aa k £ p+ iq kv + r J 

which, separated in real and imaginary parts, gives 
the quantities R aa and I aa . Similarly, the remaining 
quantities R and I are obtained. They are all func- 
tions of k or C{k). The terms with bars R aa , R ta , 

and R ch are seen to be the only ones containing tin; 
unknown X. Tin; quantities V. and A will be defined 
shortly. The quantities R find 1 are given in the 
following list: 
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'A.— ^. 1+ iur+i2(l+«)[(i-a)e-^] (1) 

1 Rt>=— A si + -£ 1 A/i + £-(a + 2)^“®“ 2 (“) 

R a /i= —Ah\ + ^2^a 4 ^<7 -( 3 ): 

(4) 

■ S. ( --B„+^B»+SVX-^[j’„G-2r,^] (5) 


• Complex cubic equation in X: 

r ^a^/3'^cft4*^^/i-4ua‘f _ ^h^a'46/s)^Y 2 

4“ v^aA/aa 4" Afft/S 4~ ®hM C h) X-\-D~Q (XXI ) 

Case 3, |3) : 

S2 a tt^X 2 “l-- (fi a ^i60-i-^0'^aa)X-]- jM ( h — 0 (XXII) 

Case 2, (fi, h)\ 

J2^aX 2 + {SlpA t ) l -\-Sl fi Ai)p)X-\~M aa =Q (XXIII) 

Case 1, {h, a ) : 


t? T> 1 -7 12 /if 

^ 6 *=~2?m jt ^ ti 

|«^-CV- 4 2[(I-a)(?-|F] 

'fl« # ==-c s ,-II[r 11 (?-2r 1 oj?] 

ftch — — (7/,1 + flfl.Y— ^2(? 

\ ft 

/„--i[2( 0 + i)!(|- 0 )F+ 4 (?}-^] 

' ^!- - ®( a + £> ( T " F+2 E Tl » e ) - A *] 
/„--I 2 ( 0 + i)F 

’W[T{G-O f+ H+*-] 

' 1 / »-5[^( r » ir + 2 l y “ e )+ B «] 

t 1 3^12 rr> 

[ I »' m tr F 

’ 7 "-f[ 2 {(!-«>+! ®} +<7 -]. 

* / -“iH!K r,,F+2 r r " fi, ) +c -] 


( 6 ) 

(7) 

( 8 ) 
(9) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 
(19) 


The solution as given by the three-row determinant 
shall be written explicitly in X. We are immediately 
able to put down for the general case a cubic equation 
in X with complex coefficients and can easily segregate 
the three subcases. The quantity D is as before the 
value of the determinant, but with the term containing 
X missing. The quantities M aa , M h $> and M eh are 
the minors of the elements in the diagonal squares 


aa , 6/3, and ch, respectively. They are expressed ex- 
plicitly in terms of R and I under the subcases treated 

in the following 

paragraphs. 

j A aa + Q a X 

A,:, A ah ! 

L)— | Af,a 

h/.H” ilgX A bh =0 

j A cat 

Acs 1 


where .4 aa[ =.Raa-H/a« etc. 


fi*fi a X 2 4' (Slj*^aa + &<^ c fc)X+iWwJ==0 






il^X — 


C a = / u a r a \ 2 1 / br T co f V 
*W<c"WJ K\vk ) 

C s { upr $ V 1 / br T 0>r V 

k 2 Mv 2 n \u T r t J K\vkJ 

Cnb* Y^Vl (br T u T V 
k 2 Mv 2 K \u,r T J ti\vk ) 


and finally 


(XXIV) 


We are at liberty to introduce the reference param- 
eters u r and r Tt and the convention adopted is: u r is 
the last u in cyclic order in each of the subcases 3, 2, 
and 1. 

Then ) and fi* +1 = l, thus for 

V^n+l^n+l/ 

Case 3, and 

Case 2, and ^/»=1 

Case l, f and Sl a — 1 
\o>«r a/ 


To treat the general case of three degrees of freedom 
(equation (XXI)), it is observed that the real part 
of the equation is of third degree while the imaginary 
part furnishes an equation of second degree. The 
problem is to find values of X satisfying both equ a- 
tions. We shall adopt the following procedure: Plot 

graphically X against ^ for both equations. The points 
of intersection are the solutions. We are only con- 
cerned with positive values of ^ and positive values of 

X . Observe that we do not have to solve for k , but 
may reverse the process by choosing a number of 
values of k and solve for X. The plotting of X 

against -g for the second-degree equation is simple 

enough, win* re as the : : • : ,o is somewhat more 

laborious .. *r the i ,- 1 iu.it ion. Ho vvever, 

the general case is of less pi cyclical importance than 
are the three subcases. ., The equation simplifies con- 
siderably, becoming of second degree in X. 
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We shall now proceed to consider these three sub- 
cases. By virtue of the cyclic arrangement, we need 
only consider the first case (<*, /S). The complex 
quadratic equations (XXII)-(XXIV) all resolve 
themselves into two independent statements, which 
we shall for convenience denote “Imaginary equa- 
tion” and “Real equation”, the former being of first 
and the latter of second degree in X. All constants 
are to be resolved into their real and imaginary parts, 
denoted by an upper index R or I, respectively. 

Let M aa —M R aa -\-iM I aa and let similar expressions 
denote M bfi and M ch 

Case 3, (a,/3). Separating equation (XXII) we obtain. 
(1) Imaginary equation: 

(fl albp T flfl/oa) X T M r c >, — 0 

y M 7 cft 

®al bfiV 

(2^ Real equation: 

tiatyX 3 -^ (i'l a Rt,0-{-^fRaa)X-}-M R ch — O 
Eliminating X we get 

Sl a &p(M T c ),) 2 — (tl a Rw-{-SlpRaa) Iaa)M*ch 

*T M R ch (flaiftlS + ^ 0 laa) 2 = 0 

By the convention adopted we have in this case: 

andc ^ 1 

Arranging the equation in powers of ft* we have: 

vi + 8 a m\ h y 

— A / 7 , ch {Raalbfr J rIa*Rbfi) J r 2 M R chIaaI 6/s] 

+ [- M' eh R aa I aa + M\ h I aa 2 } = 0 

But we have 

(AfT*) 4 — M. 1 ch(RaaIb$~ i rI<iJRbei) 

^1 Ch\.RaaIbfl Ra$I ba I Rb@Iaa Rbalafi Raal bff' ~Rb$I m>] 

= A/ 7 ch(Rat)Iba T" I aflRba) 

Finally, the equation for Case 3 (a, ff) becomes: 

Sla(M R chl b $ 2 — M 1 ckR bfil b fi) -f-ft a [ — M 1 ch{R a$I ba~\~ I a$Rba) 

+ 2M R eh I aa I, 6(3 ] + M R ch I aa 2 — M^^Raalaa = 0 (XXV) 
where 


Case 2, (0, h) fl„=(^) rj ! S!»= 1 

+2M*JvI„]+MW-MiJl M! I ll ^0 ( XXVI) 
where MV = R tf R, h — R h «R ci - / M / t ,+ 

A i„a Rb&I ch RbJcd “F~ I bfllch I bh Red 

Case 1, (A, a) ai f =« a ft/,= f— ^ — 2 ft„=l 

\i» a / r a 

ttliMWaJ A ibpRaalaa) T &h{~~ Mbji(R cal ah~T I Ca R ah) 
+2M% & I c J aa \+M%IJ~M b fi ch I ( ^Q (XXVII) 
where M% ^R ch R aa — R ca R ah — I ch I aa 4- I ca I afl 

A I b 0 — Rchl act R cal ah -f* I chRaa IcaRah 

Equations (XXV), (XXVI), and (XXVII) thus 
give the solutions of the cases: torsion-aileron , aileron- 
deflection, and deflection-torsion, respectively. The 
quantity ft may immediately be plotted against 

^ for any value of the independent parameters. 

The coefficients in the equations are all given in terms 
of R and I, which quantities have been defined above. 
Routine calculations and graphs giving ft against 

^ are contained in Appendix I and Appendix II. 

Knowing related values of ft and p X is immediately 

expressed as a function of ft by means of the first- 
degree equation. The definition of X and ft for each 
subcase is given above. The cyclic arrangement of 
all quantities is very convenient as it permits identical 
treatment of the three subcases. 

It shall finally be repeated that the above solutions 
represent the border case of unstable equilibrium. 
The plot of X against ft gives a boundary curve between 
the stable and the unstable regions in the Aft plane. 

It is preferable, however, to plot the quantity p 

instead of X, since this quantity is proportional to the 
square of the flutter speed. The .stable area can easily 
be identified by inspection as it will contain the axis 

p ^ ^ ie combination is stable for zero velocity. 


AI ch RaaRbff Ra&Rba Iaalb0~ f” lafilba 

Af ch Raalbff Ra&I ba~V IgaR b& Ia$Rba _ _ , 

Langley Memorial Aeronautical Laboratory, 

The remaining cases may be obtained by cyclic National Advisory Committee for Aeronautics, 
rearrangement: ! Langley Field, Ya., May 2, 1934. 



APPENDIX I 

PROCEDURE IN SOLVING NUMERICAL EXAMPLES 


(1) Determine the R’s and J’s , nine of each for a 
major case of three degrees of freedom, or those per- 
taining to a particular subcase, 4 R’s and 4 I’s. Refer 
to the following for the R’s and I’s involved in each 
caser ' • - \ 

The numerals 1 to 9 and 11 to 19 are used for con- 
venience. 

(Major case) Three 
degrees of freedom 


1 

Raa 

7 aa 

11 

2' 

RaP 

lap 

12 

3 

Rah 

I ah 

13 . 

4 

Rf>a 

7 ba 

14 

5 

Rb& 

I bp 

15 

6 

Rbh 

I bh 

16 

7 

Rea 

lea 

17 

8 

B c & 

I* 

18 

9 

Be h . 

I ch 

19 


(Case 3) Torsional- 
aileron (a, /3) 

1 

Raa 

I aa 

11 

2 

Rap 

lap 

12 

4 

Rba 

I ba 

14 

5 

Rb$ 

IbP 

15 


(Case 2) Aileron- 
deflection (0, h) 

5 

RbP 

IbP 

15 

6 

Rbh 

I bh 

16 

8 

Re p 

IcP 

18 

9 

Rch 

I ch 

19 


(Case 1) Deflectiom 
torsion ( h , a) 

7 

Rea 

I ca 

17 

9 

Bch 

Ich 

19 

1 

Raa 

I aa 

11 

3 

Rah 

I ah 

13 


tions of the two independent parameters a and c only. 6 
The formulas are given in the following list. 


2 1 


( 1 ) 

. _11J 
R a$ ~k ttI 

[(r.+r»)i+(<.+|)(rnfl-|r^){ 

(2) 

B"«» = j2| 


(3) 


(4) 

Dll 1 

K bfi ~ ~~ Tc 


(5) 

R" — n 

R bh ~ 

(6) 

Till 1 

K ca~ k 


> (7) 


(8) 

Tilt 1 

R eh~ k 

2 G ’ 

■ (9) : 

ha- ~2 

( a+ I)|(l _a )*’ + r Gf } + ^~ a 

(11) 

u=-\ 

{(a+|)(r„r+lii„(?)+2 2 , 

(12) 

< 



I ah— 2 

(," ! l) F 

(13) 

it i 

[(i_ a ) / . + i G l + I( p -r 1 -ir ) ) 

(14) 

Where p = 



J __L 

lb& 2tt 2 


(15) 


F 

(16) 

L. = 2j( 

|-«)r+|6 ! }+ 1 

(17) 


7’„F+|r,off)-r.} 

(18) 


It has been found convenient to split the h s in two j 
parts R=R'+R", the former being independent of j 

the argument k . The quantities I and R are func- j 


T —oF 

1 ch - 1 


( 19 ) 


s The quantities / given iu tbs appendix and used in the mowing caieu 

are seen to differ from the I’s given in the body of the paper by the factor 

may be noticed that this factor drops out in uhe first-degree equations. 


k 


iuns 


It 
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: • ! Case 3 


Choosing certain values of a and c and employing 
the values of the T’s given by the formulas of the report 
(p. 5) or in table I and also using the values of F and 
G (formulas (XII) and (XIII)) or table II, we evaluate 

the quantities / and R" for a certain nu\nber of ^ 

values. The results of this evaluation are given in 
tables III and IV, which have been worked out for 
a = 0,— 0.2, and— 0.4, and for c = 0.5 and 0. The 

range of ^is from 0 to 40. These tables save the work 

of calculating the I’s and R" ’s for almost all cases of 
practical importance. Interpolation may be used for 
intermediate values. This leaves the quantities R' to 

be determined. These, being independent of p are as 

a result easy to obtain. Their values, using the same 
system of numbers for identification, and referring to 
the definition of the original independent variables on 
pages 9 and 10, are given as follows: 


R’r 


J?' ^8 


7-(g +a2 ) 

(c-a) 


^ + ^+(c-a)^ 

K 7T 7T 


( 1 ) 

(2) 


K 

(3) 

jR' 6 a = saine as R ' a 0 

(4) 

7 ? r ^*0 i I 7 1 

(5) 

K 7T 

Si 

+ 

Si * 
1 

II 

■c 

ft? 

( 6 ) 

j£' ca =same as R' ah 

(7) 

R' cl 9 = same as R\ h 

( 8 ) 

1 

K 

(9) 


Because of the symmetrical arrangement in the 
determinant, the 9 quantities are seen to reduce to 
6 quantities to be calculated. It is very fortunate, 
indeed, that all the remaining variables segregate them- 
selves in the 6 values of R' which are independent of p 

while the more complicated I and R" are functions 
solely of c and a. In order to solve any problem it is 
therefore only necessary to refer to tables III and IV 
and then to calculate the 6 values of R' . 

The quantities (1) to (9) and (11) to (19) thus 

having been determined, the plot of 12 against p which 

constitutes our method of solution, is obtained by 
solving the equation -{-c-=0. The constants. 

a, b, and c are obtained automatically by computation 
according to the following scheme: 


Find products 1.5 


2.4 
1 


Y -V 

11.15 


15 


12.14 


men hF^=1.5— 2.4-p(ll. 15-12.14) C - 

Find products 1.15 2.14 11.5 

Then M r ch =l. 15-2. 14 + 11. 5-12.4 
and a=M R ch ( 1 5) 2 — AF eh (5.15) 

5 = ~M T ch (2 . 1 4 -f- 1 2 .4) + 2M* C ,, (1 1.15) 
c—M R ch (l l) 2 — AFcaQ.II) Find Sl a 

Solution:* — ^\ 5 J± U 
A M ! ch 


12.4 


Similarly 


Case 2 



5.9 6.8 

15.19 

16.18 

M\„ 

-5.9—6.8 — 4(15.19 

-16.18) 



5.19 6.18 

15.9 

16.8 


-5.19-6.184-15.9- 

16.8 



a = M R aa ( 1 9) 2 — M f aa (9.19) 
b——M I aa (6.18+16.8) 

+ 2M* aa (6.18-f-16.8) 
c^M R aa (l5) 2 -M I aa (5.15) 
1 12^(19)4*15 

X M I aa 


Find 0-3 


and 


Case 1 


9.1 


M R 


b8 z 


7.3 

19.11 

17.13 

-p(19.11 

-17.13) 


7.13 

19.1 

17.3 


9.11 

M T b 0~9.11 —7.13 + 19.1-- 17. 3 

b= -M 1 ^. 13 + 17.3) + 2M R ^(19. 11) 
c —M R b& {\9)— M / {, /3 (9.19) Find £2, 
1 12,(11)4-19 




Z 

M\t 




S2« 

is 

defined as ( 

( u a r a y 

for 

case 

3; 








£2^ 

is 

defined as I 

ft?)’ 

for 

case 

2; and 

£2/» 

is 

defined as | 

fey 

for 

case 

1. 


1 . / vk V 

The quantity y ,"T ) by definition. 

Since both 12 and are calculated for each value of 

p we may plot p directly as a function of 12. This 

quantity, which is proportional to the square of the 


T 


quantity, viz, j y A 


2 —> and will 
bo} r r r 


toot of the above 

denote this 
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quantity by F, which we shall term the “flutter factor/' 
The flutter velocity is consequently obtained as 


v~F 


boir^r 

w 


Since F is nondimensional, the quantity 


bo) r T f 

V* 


must 


obviously be a velocity. It is useful to establish the 
significance of this velocity, with reference to which 
the flutter speed, so to speak, is measured. Observing 
irob 2 

■■ and that the stiffness in case 1 is given by 


that 


v 


C a 


Mb 2 r c 


this reference velocity may be written: 
bo3 a r a l /O, 


VjT 


V* 


"W 




or 


rpv R 2 b 2 =C a 


The velocity v R is thus the velocity at which the total 

force on the airfoil irpv R 2b attacking with an arm | 

equals the torsional stiffness C a of the wing. This 
statement means, in case 1, that the reference velocity 
used is equal to the “divergence" velocity obtained 
with the torsional axis in the middle of the chord. This 
velocity is considerably smaller than the usual diver- 
gence velocity, which may be expressed as 

1 

Vd = Vr ~ 

2 + a 

where a ranges from 0 to — We may thus express 
the flutter velocity as 

v f ~v h F 


In case 3 the reference velocity has a similar signifi- 
cance, that is, it i.s the velocity at which the entire lift of 

the airfoil attacking with a leverage 7 / b equals numeri- 
cally the torsional stiffness Cp of the aileron or movable 
tail surface. 

In case 2, no suitable or useful significance of the 
reference velocity is available. 

TABLE I.— VALUES OF T 




c='A 

e=*0 

c = -K 

e= — 1 

Tt 

0 

-0. 1259 

-0. 6067 

-1. 6967 

-3.1416 

Tj 

0 

-0. 2103 

-1.5707 

-4. 8356 

-9.8697 

T} , 

0 

05313 

-.8084 

-3. 8375 

-11. 1034 

T t 

0 

~. 6142 

-1.5708 

-1. 6614 

-3. 1416 

TV 

c 

-.9398 

-3. 4674 

—6. 9503 

-9. 8697 

Tt 

0 

-0. 2103 

-1.5707 

-4. 8356 

-9.8697 

T? 

0 

.0132 

-. 1964 

-1. 1913 

-3. 5343 

Tt 

0 

.0903 

-.3333 

-1. 4805 

-3. 1416 

Tit - 

0 

1. 9132 

2. 5708 

2. 9604 

3. 1416 

T11 

0 

1.2990 

3. 5708 

6. 3538 

9. 4248 

Tn 

0 

. 07066 

.4292 

1. 2990 

3. 1416 


TABLE II.— TABLE OF THE BESSEL FUNCTIONS / 1 
Yo, Yi AND THE FUNCTIONS F AND G 

..... ■/i( J r l +yo)-fVi(Vi-./o) 


-0(A) 


Y,(Ji + Vc)- J,(Y,- Jo) 
(Ui+ Vo) J *M Yi~ Uo) J 


k 

1 

A 

Jo 

J\ 

Vo 

Yi 

F 

— G 

<0 

0 





0. 5000 

0 

10 

Ho 

-0.2459 

0. 0435 

0, 0557 

0.2490 

. 5006 

0. 0126 

6 

H 

. 1506 

-. 2767 

-.2882 

-. 1750 

. 5018 

.0207 

4 

H 

-. 3972 

-.0660 

-.0170 

.3979 

. 5037 

.0305 

2 

h 

.2239 

. 5767 

.5104 

-. 1071 

. 5129 

.0577 

1 

1 

.7652 

. 4401 

.0882 

~. 7813 

.5395 

. 1003 

.8 

1H 

.8463 

.3688 

-.0868 

-. 9780 

. 5541 

. 1165 

.6 

m 

.9120 

.2867 

-.3085 

-1.2604 

.5788 

. 1378 

.5 

2 

.9385 

.2423 

-. 4444 

- 1. 4714 

.6030 

. 151 

.4 

2 H 

. 9604 

. 1960 

-. 6060 

- 1 . 7808 

. 6245 

. 166 

.3 

3H 

.9776 

. 1483 

-. 8072 

— 2 2929 

.66.50 

. 180 

.2 

5 

.9900 

.0995 

-1. 0810 

—3. 3235 

. 7276 

. 1886 

.1 

10 

. 9975 

.0499 

-1.5342 

-7.0317 

.8457 

. 1626 

. 05 

20 





. 911 

. 132 

. 025 

40 





.965 

1.000 

. 090 

0 

ao 





0 
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TABLE III.— VALUES OF R 


1 

k 

0 

Ho 

H 

M 

H 

1 

in 

1H 

\7 

2 

m 

3>$ 

5 

10 

20 

40 


e 

a 







0 

0 

-0.00664 

■HI 

BBS 

ffiM 

-0. 58965 

■HI 

-1.72330 

—2 56300 

mm, 

-7. 68720 




wmm 

R".m 

« 

-.2 

0 


ioIuIjaI 

-. 02208 


-. 36586 

-.58061 


-1.57400 

-2.51580 


-11. 31010 

-51. 42490 





-.4 

0 

sSuia 

BLucl 


BLaii 

-. 12595 

-. 19936 


~. 53676 

-.85520 

BU 


-17. 20670 

-73. 35620 

-305. 9280 


IP 


0 




P*iftrj! 

iMKRiffi 

-.29384 

-.56223 

-.87212 

-1. 43983 

-2.84988 


-38. 296.50 

-172.36360 

-741. 7972 



-.2 

0 

iHiffrSTi; 


tiHvijif:* 

BKuSw 

HfSttfyS 

. 02266 

. 01629 

— 


-.29517 

-1. 29480 


-52. 49020 

-241. 3664 

ft”.* 


-.4 

0 

.00222 

.00617 

KMi 

PEa 

.21861 

.33914 

. 59199 

.84414 


2. 25914 

4. 87340 

17.80470 

67. 38320 

259. 0648 


0 

0 

1 MfeWgSs! 



iBfffR? 


.08876 

. 12176 

.12260 

MR3?R 


-.93535 

-10. 48970 

-59. 16180 

-268. 7236 


as 

-.2 

0 

.00214 

Mjhj- 


.05278 


. 31065 


.73222 

1, 10233 

1.81135 

3. 55230 

10. 14740 

31.49620 

» i MT il 



-.4 

0 

.00347 




. 34361 

. 53463 

. 94336 

1. 34762 


3. 66913 

8. 08235 

30. 97980 

120. 69760 

475. 2592 

R".k 

(*) 

0 

0 


-.00345 


ISSEUE 

PJHWIfj 


-. 22470 


MR 


IBSt'lkSS 

-1. 62600 

HUB BWj 


-.2 

0 


-.00207 

^SImsKS 

-.01734 



13482 

-. 18120 

Bill'll' l 

^^SRCtrifr 

-.56580 

-.97560 

(it 



WM 

ES 

0 


—.00334 

Bum 


as 




HKieiiuii 

-.06683 

-. 10030 

-.20060 

-.40120 

-.8024 


■P 


0 


.00214 


. 01949 

. 08055 

.12821 

.23541 

.35010 

.56143 

1.05008 

2.54920 

11. 66330 


208.7520 


||K 


0 

MBtuIiH*! 


SHliiticiX 


. 08329 

. 13219 

. 24169 

.35836 

. 57276 

1. 06650 

2.57430 



208. 8500 



EL 

0 

BaUfiit 


mm 

.02106 


. 13616 

.24796 

.36661 

maim 

1.08286 

2. 60069 

11. 75220 


208. 9490 




0 

HRf 


, ^ *1 r jSKS 

. 00321 

.01327 

.02112 


. 057S7 

.09248 

. 17296 

.41988 

1.92110 




0 . s 

Hr 

0 

dBvijMk 




.01372 

.02177 



.09434 

. 17566 

. 42413 

1.92840 

;^Klw 




BL 

0 

sK’to.ot 1 


SSS 


.01417 

.02243 



. 09621 

. 17836 

.42837 

1. 93575 

8. 25246 

34, 4169 

R"bfi 

a 


0 

.00124 


. 00772 

. 03101 


. 19830 

.35807 

ilHKfeCT- 

■n 


3. 54970 

15. 35120 

64. 02240 

263. 2340 


(») 

0| 

.00031 



.00785 


.04680 



aBESs B 


. 84970 

3. 55050 

14. 56740 

59. 3188 

RTkk 

c 


0 

.00017 

1H7S| 


. 00394 


.01989 



MB 


. 12881 

.22211 

.36062 

,4918 

.6 

( j ) 

0 

IBElUiis] 




. 00226 

,00328 



.00934 


. 02122 

• 03650 

. 05940 

.0810 


« 

K 

0 

.01128 

.03132 


spt.S 

1. 17930 


3. 44670 

mbw&v 


15. 37450 



mcnxiss 

3, 056. .4000 

R ,e 

E 

0 

.01178 



. 29084 

1.21954 

1. 93540 

6. r -i860 



15. 61440 


17). 41640 


3, 057. 8400 



iBSu 

0 

.01228 



. 39838 

1. 25950 


3. e.;. yj 

5. 307 60 

8. 55200 


38. 07740 


|^(||2gas2 

3, 059. 2800 

RT'fi 

HP 

<*) 

O' 




-. 24'266 

-1.0056) 

-1. £8246 

-2.S9371 

m 


-15.49965 

NiKKcliMl 

- 140. 26370 


-2,502.3470 

.6 


0 




. 16810 

IHLiMjJ 

1.11453 

2. 05320 


9. 24438 

KgjymjJ 


|Bt2£2S 

1,881.9900 

P" t k 


(*) 


.00250 

.00690 

.01420 

. 05780 

.20060 

.29120 

.44940 


.83000 

1.20000 

1, 88600 


5.28000 

7.2000 


1 Independent ot c. * Independent of a. 

TABLE IV.— VALUES OF I 



■ 

Ho 


M 

H 

1 

W 

m 


2*S 

3H 

5 

10 

20 

40 

■ 

a 

a 

■ 


la. 

o) 

0 

0.25000 

0. 25096 

0. 25255 

0. 25578 

0.27240 





0.60275 

0.78750 

1.07920 

1. 70320 

2.68450 

3.61750 

-.2 

.49000 

.49050 

.49131 

. 49302 

.50189 

. 53359 

.55464 

. 59472 

.62794 

.68671 

.78070 

.96021 

1. 32940 

1. 90140 

2. 45470 



-.4 

.81000 

.81014 

.81037 

. 81086 

.81145 

.82395 

.82938 

.84176 

.85186 

. 87059 

.00030 

. 95763 

1.07300 

1.26400 

1.44630 


Wm 

0 

. 17805 

. 17874 

. 17985 

. 18219 

. 19433 

.23768 

.26645 

. 32132 

.36664 


. 57526 

. 82035 

1.31213 

2. 10476 

2.85963 


■ 

-.2 

.39170 

. 39212 

.39278 

.39418 

. 40147 

.42748 

. 44474 

.47761 

.50485 



.77708 

1.07215 

1.54773 

2. 00065 

I.fi 


-.4 

.60531 

.60545 

,60567 

.60614 

.60857 

. 61724 

.62299 

.63395 

.64303 

.65908 

.68475 

. 73377 

.82313 

.99065 

1.14163 


0 

. 13252 

. 13317 

. 13425 

. 13640 

. 14742 

. 18544 

.20914 

.25611 

. 29514 

. 35951 

.46379 

. 65973 

1.05124 

1. 65524 

2.22869 


-.2 

.21297 

.21336 

. 21401 

.21530 

. 22191 

.24472 

.25894 

.28712 

.31054 

.34916 

.41173 

.52929 

.76420 

1. 12651 

1. 47067 


a 

-.4 

.29342 

.29354 

.29376 

.29419 

.29640 

.30400 

. 30891 

.31813 

. 32594 

.33881 

. 35966 

.39884 

.47714 

. 59792 

.71260 

I.k 

o) 

Kfl 

-.50000 

-.50060 



-.51290 


-.55410 

-. 57880 

w]B 

62450 

-.66500 

-.72760 

-. 84570 

-.94100 

—.96500 

31 

-.30000 

-.30036 


- 30222 

-. 30774 

sHMaftl 

-.33246 

-. 34728 


37470 

-. 39900 

-. 43656 

-. 50752 

-.56460 

-.57900 



n 

10000 

-.10012 



-. 10268 

-. 10790 

11082 

-. 11576 

12060 


13300 

-. 14552 

-.16914 

-. 18220 

-. 19300 



0 

.39023 

.39010 

.38988 

. 38944 

.38717 

.37923 

. 37404 

.36424 


.34204 

.31954 

.27696 

. 19172 

.05766 

— . 06080 


0 

-.2 

.40389 

.40378 

.40359 

.40320 

.40119 

. 39397 

. 38918 

.38005 

. 37249 

.35911 

. 33771 

.29683 

. 21469 

.08255 

-.04344 



-.4 

. 41755 

. 41746 

. 41730 

. 41697 

.41520 

. 40871 

.40432 

.39586 

.38896 

.37617 

. 35598 

. 31671 

.23793 

.10744 

-. 01707 


m 

o 

.07438 

.07435 

.07433 

.07424 

.07387 

.07256 

.07171 

.07009 

.06874 

.06644 

.06273 

. 05572 

.04168 

.01960 

-.00327 


-.2 

101 


. 07658 

. 07651 

.07618 

.07469 

.07420 

.07270 

. 07145 

.06925 

.06572 

.05899 

. 04518 

.02370 

.00295 


a 

-.4 

. 07887 


.07882 

.07867 

. 07848 

.07741 

.07668 

.07529 

.07416 

.07205 

.06871 

.06226 

.04928 

.02779 

.00728 

fbfi 

m 

( J ) 

.32297 

.32288 

R&31 

.32241 

.32075 

.31483 

. 31090 


.29721 

.28625 

.26872 

.23524 

.16806 

. 05979 

04333 

■9 


.04270 

.04270 

.04270 


.04240 

.04150 

.04095 


.03904 

.03760 

.03386 

.03080 

.02200 

.00845 

-.00470 

hk 

0 

(*) 

.06830 

.06840 

.06850 



Miiilf?!! 

. 07570 

. 07910 

.08240 

.08530 

.09080 

.09940 

.11550 

.12440 

. 13180 

.5 


.01125 

. 0L12C 

.01129 



.01214 

. 01247 

. 01302 

.01357 

. 01405 

.01496 

.01637 

.01903 

.02117 

.02171 

■■j 

0) 




1. 49490 

1. 48844 

1.45520 



1. 11940 


WEEtSl 







-.2 


1.69832 

1. 69562 

1. 68992 

.1,66036 

1. 55*170 

1.48454 

1.35092 


1.04430 


. 13264 

-1.06802 


-4. 84900 

■ 


-.4 


1.89856 

1. 89634 


1. 88552 


1.70618 

1.58244 


1.31410 



-.72974 



HB 

0 

(*) 



1. 06470 

1. 06000 


mjgffini 

. 891.50 

.78190 

.69110 

. 52840 

.27380 

-.21640 

-1.20010 

-2. 78550 

-4. 29530 

mm 

.5 





.39450 


.29640 

. 24640 

. 15510 



—.26030 

-. 65220 

-1.43520 

-2. 64380 



0) 

(») 



3 


1.02580 

1.07900 

1. 10820 

1. 15760 

1.20600 

1.24900 

1.33000 

1.45520 

1. 69140 

1.82200 

1. 93000 


1 Independent of e. > Independent of a. 





















































































































































































APPENDIX li 

NUMERICAL CALCULATIONS 


A number of routine examples have been worked out 
to illustrate typical results. A “standard” case has 
been chosen, represented by the following constants: 





k~0. 1 , c«= 0 . 5 , ««=— 0 . 4 , z a = 0 . 2 , 

Ta 2 ~0.25 t zg b gQ> Tp 

« a , w* variable. 


We will show the results of a numerical computation 
of the three possible subcases in succession. 


/GO 

mo 

n. 

80 

40 
0 

Figure 5.— Case 3, Torsion-aileron (a, 0): Standard case. Showfng against 



Case 3, Torsion-aileron (a,i 3): Figure 5 shows the ft a 
against ^ relation and figure 6 the final curve 

M^J against fi “ = (S^) s=40 (S) 



Figure 6. — Case 3, Torsion-aileron (or, 0): Standard case. Showing flutter factor 
F against 12„. 


Case 2 , Aileron-flexure ( 0 , h ): Figure 7 shows the 

1 . . . . / r \ 7 

% against ^ relation* and figure 8 the final curve f 


against ft 0 = 


( 


'wA : 

\ u h , 

1 160V 



• It is realized that considerable care must be exercised to get these curves reason- 
ably accurate. 


The heavy line shows the standard case, while the 
remaining curves show the effect of a change in the 

value of xg to jx and —U* 


Case 1 , Flexure-torsion ( h , a): Figure 9 shows again 



Figure 7.— Case 2, Aileron-deflection (0, ft): (a) Standard case, (b), (c), (d) indicate 
dependency oaxfi. Case (d),.Z 0 = — 0.004, reduces to a point. 

the ft/, against ^ relation and figure 10 the final result 



against «. = (~) ““(S) 


Case 1, which is of importance in the propeller theory, 
has been treated in more detail. The quantity F shown 


v 


in the figures is V* ~~V~h 

U/ a f a ” 


Figure 11 shows the dependency oil ^ 1 = 

figure 12 shows the dependency on the location of the 
axis a; figure 13 shows the dependency on the radius of 
gyration r a =r; and figure 14 shows the dependency 
on the location of the center of gravity x, for three 
different combinations of constants. 

; -. ... .■ i / - .. . i •. L I S 

Detailed discussion of the experimental work will not 
be given in this paper, but shall be reserved for a later 
report. The experiments given in the following are 
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restricted to wings of a large aspect ratio, arranged with 
two or three degrees of freedom in accordance with the 



Figure 8. — Cn,w 2, Aileron-defleetion (0, A): Final curves giving flutter factor F 
against fa corresponding to cases shown in figure 7. 


theoretical cases. The wing is free to move parallel to 
itself in a vertical direction ( h ); is equipped with an 



axis in roller bearings at (a) (fig. 2) for torsion, and 
with an aileron hinged at (c). Variable or exchange- 


able Springs restrain the wing to its equilibrium 
position. 





Figure 10.— Case !, Flexure-torsion (ft, a): Standard case. Showing llutter factor 
F against fa. 

We shall present results obtained on two wings, both 
of symmetrical cross section 12 percent thick, and with 
chord 26 = 12.7 cm, tested at 0°. 



CjJj/CJs. 


Figure II.— Case 1, Flexure-torsion (ft, «): Showing dependency of Fon^* The 
up|ver curve isexpcrimental. Airfoil withr= 1- a = — 0.4;r=0.2; 4*=.0I ; ~ variable. 

Wing A, aluminum. vbh f 1 r 11 x vimr constants: 

k — (<--■- ./■„ !•. - 1 , i"! i n, and 0.038, 

4 lb ' * ’ 

respectively; 


r a 2 =0.33 and a>„=7X2ir 
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Wing B, wood, with flap, and the constants: 

c=0.5, a=— 0.4, r a = 0.192, r a 2 =0.178, 

3^=0.019, V=0.0079, and u a kept constant 

= 17.6X2* 

The results for wing A, case 1, are given in figure 15; 
and those for wing B, cases 2 and 3, are given in figures 
16 and 17, respectively. The abscissas are the fre- 
quency ratios and the ordinates are the velocities in 
cm /sec. Compared with the theoretical results calcu- 
lated for the three test cases, there i3 an almost perfect 



Fiqcee 12.— Case 1, Flexure-torsion (h, a): Showing dependency of F on location 

„ — . . h; a variable. 

2 4 0 


of axis of rotation a. Airfoil with r- A ; x=0.2; * —A ; — « 


agreement in case 1 (fig. 15). Not only is the minimum 
velocity found near the same frequency ratio, but the 
experimental and theoretical values are, furthermore, 
very nearly alike. Very important is also the fact that 
the peculiar shape of the response curve in case 2, pre- 
dicted by the theory, repeats itself experimentally. 
The theory predicts a range of instabilities extending 
from a small value of the velocity to a definite upper 
limit. It was very gratifying to observe that the upper 
branch of the curve not only existed but that it was 
remarks. hi y u ’dilute. A small increase in speed near 
this upper limit would suffice to change the condition 
from violent flutter to complete rest, no range of trans i- 
tion being observed. The experimental cases 2 and 3 
are compared with theoretical results given by the 
dotted lines in both figures (figs. 18 and 17). 


The conclusion from the experiments is briefly that 
the general shapes of the predicted response curves re- 



o .s /. Q t.S 


r 


Figuhe 13. — Case 1, Flexure-torsion ( h , a): Showing dependency of F on the radius 
of gyration r«=r. 

A, airfoil with a= —0.4; *—A: r=0.2! —-4-’ r variable. 

4 on o 

B, airfoil with o= —0.4; «— hi r=0.2! ——1.00; r variable. 

4 on 



O .2 -4 .6 


Figub* 14.— Case 1, Flexure-torsion (h, a): Showing dependency of F on x», the 
location of the center of gravity. 

A, airfoil with r— A< a = — 0.4; k— -A; — — -4> x variable. 

I 400 W2 o 


8, a'rfail with r a — - 0.4; 

C, airfoil with r^A; a = —0.4; 


A : 1 i variable, 

4 ti- l 2 t> 

k— -A; — — i; x variablo. 


peat themselves satisfactorily. Next, that the influ- 
ence of the internal friction 7 obviously is quite appreci- 


'i This mattor is the subject of a paper now in preparation. 
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able in case 3. This could have been expected since 
the predicted velocities and thus also the air forces on 
the ailbron are very low, and no steps were taken to 
eliminate the friction in the hinge. The outline of the 
stable region is rather vague, and the wing is subject 



Fjoukb 15. — Case 1. Wing A. Theoretical and experimental curves giving flutter 
velocity p against frequency ratio — ■ Deflection-torsion. 

OS a 

to temporary vibrations at much lower speeds than 
that at which the violent flutter starts. The above 
experiments are seen to refer to cases of exaggerated 
unbalance, and therefore of low flutter speeds. It is 
evident that the internal friction is less important at 
larger velocities. The friction does in all cases increase 
the speed at which flutter starts. 



0 .2 .4 .6 .6 LO 1.2 t.4 LG LQ 

Figure 16.— -Case 2. Wing B. Theoretical and experimental curves giving flutter 
velocity v against frequency ratio — - Aileron-deflection (£, h). 



0 0.4 0.8 L2 LG 2.0 2.4 2.8 


Wc/iOfl 

Figure 17. — Case 3. Theoretical curve giving flutter velocity against the fre- 
quency ratio The experimental unstable area is indefinite due to the im- 

OSff 

portance of internal friction at very small velocities. Torsion-aileron (<*, ff). 
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